/ Celestin Apprentice 4 / Apprentice-Release4.iso / Languages / RLaB 1.18c / manual / rank.tex

< prev

next >

MacOS 8.1

Win98

DOS

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 5c 73 65 63 74 69 6f 6e | 7b 72 61 6e 6b 20 2d 2d |\section|{rank --|
|00000010| 2d 20 52 61 6e 6b 20 6f | 66 20 61 20 4d 61 74 72 |- Rank o|f a Matr|
|00000020| 69 78 7d 0d 5c 6c 61 62 | 65 6c 7b 66 6e 3a 72 61 |ix}.\lab|el{fn:ra|
|00000030| 6e 6b 7d 0d 5c 69 6e 64 | 65 78 7b 72 61 6e 6b 28 |nk}.\ind|ex{rank(|
|00000040| 29 7d 0d 5c 69 6e 64 65 | 78 7b 72 61 6e 6b 20 6f |)}.\inde|x{rank o|
|00000050| 66 20 61 20 6d 61 74 72 | 69 78 7d 0d 5c 69 6e 64 |f a matr|ix}.\ind|
|00000060| 65 78 7b 6d 61 74 72 69 | 78 21 72 61 6e 6b 7d 0d |ex{matri|x!rank}.|
|00000070| 5c 62 65 67 69 6e 7b 72 | 61 69 6c 7d 0d 52 61 6e |\begin{r|ail}.Ran|
|00000080| 6b 20 3a 20 27 72 61 6e | 6b 27 20 27 28 27 20 4d |k : 'ran|k' '(' M|
|00000090| 61 74 72 69 78 20 28 27 | 2c 27 20 53 63 61 6c 61 |atrix ('|,' Scala|
|000000a0| 72 20 29 20 3f 20 27 29 | 27 20 3b 0d 5c 65 6e 64 |r ) ? ')|' ;.\end|
|000000b0| 7b 72 61 69 6c 7d 0d 0d | 5c 76 65 72 62 2b 72 61 |{rail}..|\verb+ra|
|000000c0| 6e 6b 2b 20 63 61 6c 63 | 75 6c 61 74 65 73 20 74 |nk+ calc|ulates t|
|000000d0| 68 65 20 72 61 6e 6b 20 | 6f 66 20 74 68 65 20 61 |he rank |of the a|
|000000e0| 72 67 75 6d 65 6e 74 20 | 6d 61 74 72 69 78 2e 20 |rgument |matrix. |
|000000f0| 20 54 68 65 20 72 61 6e | 6b 0d 6f 66 20 61 20 6d | The ran|k.of a m|
|00000100| 61 74 72 69 78 20 69 73 | 20 74 68 65 20 6e 75 6d |atrix is| the num|
|00000110| 62 65 72 20 6f 66 20 6e | 6f 6e 2d 7a 65 72 6f 20 |ber of n|on-zero |
|00000120| 73 69 6e 67 75 6c 61 72 | 20 76 61 6c 75 65 73 20 |singular| values |
|00000130| 6f 66 20 74 68 61 74 20 | 6d 61 74 72 69 78 2e 0d |of that |matrix..|
|00000140| 48 6f 77 65 76 65 72 20 | 73 69 6e 63 65 20 6e 75 |However |since nu|
|00000150| 6d 65 72 69 63 20 63 6f | 6d 70 75 74 61 74 69 6f |meric co|mputatio|
|00000160| 6e 20 69 6e 74 72 6f 64 | 75 63 65 73 20 65 72 72 |n introd|uces err|
|00000170| 6f 72 73 2c 20 74 68 65 | 20 73 69 6e 67 75 6c 61 |ors, the| singula|
|00000180| 72 0d 76 61 6c 75 65 73 | 20 61 72 65 20 61 63 74 |r.values| are act|
|00000190| 75 61 6c 6c 79 20 74 65 | 73 74 65 64 20 61 67 61 |ually te|sted aga|
|000001a0| 69 6e 73 74 20 61 20 73 | 65 74 20 74 6f 6c 65 72 |inst a s|et toler|
|000001b0| 61 6e 63 65 20 76 61 6c | 75 65 2e 20 20 54 68 65 |ance val|ue.  The|
|000001c0| 0d 6f 70 74 69 6f 6e 61 | 6c 20 73 63 61 6c 61 72 |.optiona|l scalar|
|000001d0| 20 61 72 67 75 6d 65 6e | 74 2c 20 69 66 20 70 72 | argumen|t, if pr|
|000001e0| 65 73 65 6e 74 2c 20 73 | 70 65 63 69 66 69 65 73 |esent, s|pecifies|
|000001f0| 20 74 68 65 20 74 6f 6c | 65 72 61 6e 63 65 2c 0d | the tol|erance,.|
|00000200| 6f 74 68 65 72 77 69 73 | 65 20 74 68 65 20 64 65 |otherwis|e the de|
|00000210| 66 61 75 6c 74 20 6f 66 | 20 74 68 65 20 74 72 69 |fault of| the tri|
|00000220| 70 6c 65 2d 70 72 6f 64 | 75 63 74 20 6f 66 20 74 |ple-prod|uct of t|
|00000230| 68 65 20 67 72 65 61 74 | 65 72 20 64 69 6d 65 6e |he great|er dimen|
|00000240| 73 69 6f 6e 20 6f 66 20 | 74 68 65 0d 61 72 67 75 |sion of |the.argu|
|00000250| 6d 65 6e 74 2c 20 74 68 | 65 20 32 2d 6e 6f 72 6d |ment, th|e 2-norm|
|00000260| 20 6f 66 20 74 68 65 20 | 6d 61 74 72 69 78 20 61 | of the |matrix a|
|00000270| 6e 64 20 6d 61 63 68 69 | 6e 65 20 65 70 73 69 6c |nd machi|ne epsil|
|00000280| 6f 6e 20 69 73 20 75 73 | 65 64 2e 0d 0d 24 5c 4c |on is us|ed...$\L|
|00000290| 6f 6e 67 72 69 67 68 74 | 61 72 72 6f 77 24 20 54 |ongright|arrow$ T|
|000002a0| 68 69 73 20 69 73 20 6e | 6f 74 20 61 6e 20 5c 52 |his is n|ot an \R|
|000002b0| 4c 61 42 5c 20 62 75 69 | 6c 74 2d 69 6e 20 66 75 |LaB\ bui|lt-in fu|
|000002c0| 6e 63 74 69 6f 6e 2e 20 | 20 54 68 69 73 20 66 75 |nction. | This fu|
|000002d0| 6e 63 74 69 6f 6e 0d 69 | 73 20 6e 6f 72 6d 61 6c |nction.i|s normal|
|000002e0| 6c 79 20 6c 6f 61 64 65 | 64 20 6f 6e 20 73 74 61 |ly loade|d on sta|
|000002f0| 72 74 2d 75 70 20 66 72 | 6f 6d 20 74 68 65 20 5c |rt-up fr|om the \|
|00000300| 76 65 72 62 2b 72 61 6e | 6b 2e 72 2b 20 66 69 6c |verb+ran|k.r+ fil|
|00000310| 65 20 69 6e 20 74 68 65 | 0d 73 74 61 6e 64 61 72 |e in the|.standar|
|00000320| 64 20 5c 76 65 72 62 2b | 72 6c 69 62 2b 20 64 69 |d \verb+|rlib+ di|
|00000330| 72 65 63 74 6f 72 79 2e | 20 55 73 65 20 6f 66 20 |rectory.| Use of |
|00000340| 74 68 65 20 5c 76 65 72 | 62 2b 2d 72 2b 20 6f 70 |the \ver|b+-r+ op|
|00000350| 74 69 6f 6e 2c 20 69 6e | 63 6f 72 72 65 63 74 6c |tion, in|correctl|
|00000360| 79 20 0d 73 65 74 74 69 | 6e 67 20 74 68 65 20 5c |y .setti|ng the \|
|00000370| 76 65 72 62 2b 52 4c 41 | 42 5f 4c 49 42 5f 44 49 |verb+RLA|B_LIB_DI|
|00000380| 52 2b 20 65 6e 76 69 72 | 6f 6e 6d 65 6e 74 61 6c |R+ envir|onmental|
|00000390| 20 76 61 72 69 61 62 6c | 65 2c 20 6f 72 20 6d 6f | variabl|e, or mo|
|000003a0| 64 69 66 79 69 6e 67 0d | 5c 76 65 72 62 2b 72 61 |difying.|\verb+ra|
|000003b0| 6e 6b 2e 72 2b 20 6d 61 | 79 20 6d 61 6b 65 20 74 |nk.r+ ma|y make t|
|000003c0| 68 69 73 20 66 75 6e 63 | 74 69 6f 6e 20 75 6e 61 |his func|tion una|
|000003d0| 76 61 69 6c 61 62 6c 65 | 2e 0d 0d 0d 0d 0d 5c 70 |vailable|......\p|
|000003e0| 61 72 61 67 72 61 70 68 | 7b 45 78 61 6d 70 6c 65 |aragraph|{Example|
|000003f0| 7d 0d 5c 69 6e 64 65 78 | 7b 65 78 61 6d 70 6c 65 |}.\index|{example|
|00000400| 20 6f 66 20 72 61 6e 6b | 28 29 7d 0d 5c 69 6e 64 | of rank|()}.\ind|
|00000410| 65 78 7b 72 61 6e 6b 28 | 29 21 65 78 61 6d 70 6c |ex{rank(|)!exampl|
|00000420| 65 7d 0d 5c 62 65 67 69 | 6e 7b 76 65 72 62 61 74 |e}.\begi|n{verbat|
|00000430| 69 6d 7d 0d 3e 20 61 20 | 3d 20 72 61 6e 64 28 35 |im}.> a |= rand(5|
|00000440| 2c 35 29 0d 20 61 20 3d | 0d 20 20 20 20 20 30 2e |,5). a =|.     0.|
|00000450| 34 34 20 20 20 20 20 20 | 30 2e 37 34 31 20 20 20 |44      |0.741   |
|00000460| 20 20 20 30 2e 38 34 32 | 20 20 20 20 20 20 30 2e |   0.842|      0.|
|00000470| 35 39 39 20 20 20 20 20 | 20 30 2e 32 33 39 20 20 |599     | 0.239  |
|00000480| 0d 20 20 20 20 30 2e 36 | 38 35 20 20 20 20 20 20 |.    0.6|85      |
|00000490| 30 2e 33 31 33 20 20 20 | 20 20 20 30 2e 36 35 32 |0.313   |   0.652|
|000004a0| 20 20 20 20 20 20 30 2e | 34 35 35 20 20 20 20 20 |      0.|455     |
|000004b0| 20 30 2e 34 36 35 20 20 | 0d 20 20 20 20 30 2e 31 | 0.465  |.    0.1|
|000004c0| 32 39 20 20 20 20 20 20 | 30 2e 37 35 39 20 20 20 |29      |0.759   |
|000004d0| 20 20 20 30 2e 33 38 36 | 20 20 20 20 20 20 30 2e |   0.386|      0.|
|000004e0| 36 35 38 20 20 20 20 20 | 20 30 2e 35 34 36 20 20 |658     | 0.546  |
|000004f0| 0d 20 20 20 20 20 30 2e | 37 39 20 20 20 20 20 20 |.     0.|79      |
|00000500| 30 2e 31 35 31 20 20 20 | 20 20 20 30 2e 35 38 33 |0.151   |   0.583|
|00000510| 20 20 20 20 20 20 30 2e | 32 33 32 20 20 20 20 20 |      0.|232     |
|00000520| 20 30 2e 38 30 31 20 20 | 0d 20 20 20 20 30 2e 38 | 0.801  |.    0.8|
|00000530| 36 38 20 20 20 20 20 20 | 30 2e 33 32 35 20 20 20 |68      |0.325   |
|00000540| 20 20 20 30 2e 32 30 35 | 20 20 20 20 20 20 30 2e |   0.205|      0.|
|00000550| 37 39 37 20 20 20 20 20 | 20 30 2e 32 34 37 20 20 |797     | 0.247  |
|00000560| 0d 3e 20 73 76 64 28 61 | 29 2e 73 69 67 6d 61 0d |.> svd(a|).sigma.|
|00000570| 20 73 69 67 6d 61 20 3d | 0d 20 20 20 20 20 20 32 | sigma =|.      2|
|00000580| 2e 36 20 20 20 20 20 20 | 30 2e 38 33 35 20 20 20 |.6      |0.835   |
|00000590| 20 20 20 30 2e 36 32 34 | 20 20 20 20 20 20 30 2e |   0.624|      0.|
|000005a0| 34 34 32 20 20 20 20 20 | 30 2e 30 34 37 35 20 20 |442     |0.0475  |
|000005b0| 0d 3e 20 72 61 6e 6b 28 | 61 29 0d 20 20 20 20 20 |.> rank(|a).     |
|000005c0| 20 20 20 35 0d 3e 20 62 | 20 3d 20 5b 61 2c 61 3b |   5.> b| = [a,a;|
|000005d0| 61 2c 61 5d 3b 0d 3e 20 | 73 76 64 28 62 29 2e 73 |a,a];.> |svd(b).s|
|000005e0| 69 67 6d 61 0d 20 73 69 | 67 6d 61 20 3d 0d 20 6d |igma. si|gma =. m|
|000005f0| 61 74 72 69 78 20 63 6f | 6c 75 6d 6e 73 20 31 20 |atrix co|lumns 1 |
|00000600| 74 68 72 75 20 36 0d 20 | 20 20 20 20 20 35 2e 32 |thru 6. |     5.2|
|00000610| 20 20 20 20 20 20 20 31 | 2e 36 37 20 20 20 20 20 |       1|.67     |
|00000620| 20 20 31 2e 32 35 20 20 | 20 20 20 20 30 2e 38 38 |  1.25  |    0.88|
|00000630| 34 20 20 20 20 20 20 30 | 2e 30 39 35 20 20 20 31 |4      0|.095   1|
|00000640| 2e 36 33 65 2d 31 36 20 | 20 0d 0d 20 6d 61 74 72 |.63e-16 | .. matr|
|00000650| 69 78 20 63 6f 6c 75 6d | 6e 73 20 37 20 74 68 72 |ix colum|ns 7 thr|
|00000660| 75 20 31 30 0d 20 36 2e | 31 35 65 2d 31 37 20 20 |u 10. 6.|15e-17  |
|00000670| 20 32 2e 39 33 65 2d 31 | 37 20 20 20 31 2e 30 34 | 2.93e-1|7   1.04|
|00000680| 65 2d 31 37 20 20 20 31 | 2e 34 35 65 2d 31 38 20 |e-17   1|.45e-18 |
|00000690| 20 0d 3e 20 72 61 6e 6b | 28 62 29 0d 20 20 20 20 | .> rank|(b).    |
|000006a0| 20 20 20 20 35 0d 5c 65 | 6e 64 7b 76 65 72 62 61 |    5.\e|nd{verba|
|000006b0| 74 69 6d 7d 0d 0d 53 65 | 65 20 61 6c 73 6f 3a 20 |tim}..Se|e also: |
|000006c0| 50 61 67 65 7e 5c 70 61 | 67 65 72 65 66 7b 66 6e |Page~\pa|geref{fn|
|000006d0| 3a 65 70 73 69 6c 6f 6e | 7d 20 5c 76 65 72 62 2b |:epsilon|} \verb+|
|000006e0| 65 70 73 69 6c 6f 6e 28 | 29 2b 2c 0d 50 61 67 65 |epsilon(|)+,.Page|
|000006f0| 7e 5c 70 61 67 65 72 65 | 66 7b 66 6e 3a 6e 6f 72 |~\pagere|f{fn:nor|
|00000700| 6d 7d 20 5c 76 65 72 62 | 2b 6e 6f 72 6d 28 29 2b |m} \verb|+norm()+|
|00000710| 2c 0d 50 61 67 65 7e 5c | 70 61 67 65 72 65 66 7b |,.Page~\|pageref{|
|00000720| 66 6e 3a 73 69 7a 65 7d | 20 5c 76 65 72 62 2b 73 |fn:size}| \verb+s|
|00000730| 69 7a 65 28 29 2b 2e 0d | 0d 5c 6e 65 77 70 61 67 |ize()+..|.\newpag|
|00000740| 65 0d 0d 0d 0d          |                         |e....   |        |
+--------+-------------------------+-------------------------+--------+--------+